Multi-mode mediated exchange coupling in cavity QED 



S. Filipp/'B M. Goppl, 1 J. M. Fink, 1 M. Baur, 1 R. Bianchctti, 1 L. Stcffen, 1 and A. Wallraff 1 

department of Physics, ETH Zurich, CH-8093 Zurich, Switzerland 

Microwave cavities with high quality factors enable coherent coupling of distant quantum systems. 
Virtual photons lead to a transverse exchange interaction between qubits, when they are non- 
resonant with the cavity but resonant with each other. We experimentally probe the inverse scaling 
of the inter-qubit coupling with the detuning from a cavity mode and its proportionality to the 
qubit-cavity interaction strength. We demonstrate that the enhanced coupling at higher frequencies 
is mediated by multiple higher-harmonic cavity modes. Moreover, in the case of resonant qubits, the 
symmetry properties of the system lead to an allowed two-photon transition to the doubly excited 
qubit state and the formation of a dark state. 
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I. INTRODUCTION 



Experiments on single photons coupled strongly to sin- 
gle (artificial) atoms [l| allow for in-depth studies of 
photon-atom interactions on a single particle level. This 
has first been demonstrated with individual atoms cou- 
pled to microwave [H, and later optical cavity fields 
P, Q . In solids, strong coupling has been achieved with 
quantum dots [g, 0] and superconducting circuits [1] . De- 
spite the diversity of physical realizations the coherent 
exchange of energy between photons and atoms can be 
described in all these systems by a generic model named 
after Jaynes and Cummings 0. 

In circuit quantum electrodynamics experiments, su- 
perconducting quantum circuits are coupled to single mi- 
crowave photons in a planar transmission line cavity ■ 
In this configuration, coupling strengths exceed decay 
rates by two orders of magnitude, and strong resonant 
coupling between a microwave cavity and a sing le d[T3- 
[l3j or multiple [TH, [l| superconducting qubits has been 
observed. In the case of finite detuning between a single 
qubit and a resonator mode, energy exchange between 
the individual systems is strongly suppressed due to en- 
ergy conservation. In this dispersive regime, a residual 
interaction mediated via virtual photons induces a finite 
Lamb and ac-Stark shift [l7| of the energy levels. For 
two qubits coupled to a common cavity field, the same 
mechanism leads to an interaction mediated by virtual 
photons [l8| as experimentally demonstrated |l9]. This 
coupling is similar to the J-coupling of interacting nuclear 
spins as observed in nuclear magnetic resonance experi- 
ments (e. g. [2(| HU). It is also a dominant interaction of 
double quantum dots, where the exchange splitting be- 
tween spin singlet and spin triplet state can be used to 
control a logical qubit state encoded in a two-electron 
spin-state [22|. In contrast to these local interactions, 
and also opposed to the direct coupling of superconduct- 
ing quantum circuits [23l - l29j . the coupling mediated via 



virtual resonator photons allows for a long-range interac- 
tion between two or more distant superconducting qubits. 
In the context of quantum information processing, it can 
be used to realize two-qubit gates [30l | with superconduct- 
ing qubits Hi, [Hii]. 

In this paper we measure the exchange coupling as a 
function of detuning of two qubits from a single or multi- 
ple resonator modes and characterize the symmetry prop- 
erties of the coupled system. In Section|TT]the inter-qubit 
coupling mechanism and its spectroscopic measurement 
is outlined. In Section IIIII the coupling near a single res- 
onator mode is analyzed. In Section IIVI higher harmonic 
modes of the transmission line resonator are included in 
the analysis. Section [V] and I VII describe the formation 
of a dark state at the avoided level crossing and the ob- 
servation of a two-photon transition from the ground to 
the doubly excited state that is allowed only at qubit 
resonance. 



II. EXCHANGE COUPLING MECHANISM 

In our experiments two superconducting qubits are 
dispersively coupled to a microwave cavity, see Fig. [1] 
The quantum circuits are realized as weakly anhar- 
monic transmon qubits [33j and the cavity is formed by 
a A/2 coplanar- waveguide resonator supporting several 
harmonic modes [34[. In the dispersive regime, the detun- 
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is larger than the coupling strength 



of both qubits (i — 1,2) to each resonator mode j. 
The relevant Hamiltonian 



Hi 



i=l,2 



(1) 
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is obtained by adiabatically eliminating the direct qubit- 
resonator interaction of the qubits for each harmonic 
mode dj in the Jaynes-Cummings Hamiltonian [30| • The 
first term denotes the qubit Hamiltonian with Lamb- 
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FIG. 1: (a) Coplanar waveguide resonator coupled via finger 
capacitors C K to input and output transmission lines. Two 
transmons are capacitively coupled to the resonator at its 
ends (C g ). Additional ac-signal lines are capacitively cou- 
pled to the qubits (not used in the experiments), (b) Opti- 
cal micrograph of a transmon qubit. (c) Schematics of the 
measurement setup. The state of the qubit is determined 
by measuring the transmission of the RF signal through the 
transmission line cavity modeled as an LCR oscillator. When 
a spectroscopy signal (Spec) is resonant with a qubit transi- 
tion, the resonance frequency of the cavity is shifted and the 
change in transmission amplitude is recorded at the analog- 
digital converter (ADC) after down-conversion with a local 
oscillator (LO) [351 ] . The qubit frequencies can be tuned in- 
dependently with superconducting coils (SC Coil 1/2). 

shifted transition frequencies uj g 1 } from the ground to the 
first excited state. Higher transmon levels do not play a 
role in our experiments and are therefore neglected. The 
second term in Eq. (TT]) describes the resonator modes 
with frequencies Uj = (j + l)wo, integer multiples of the 
fundamental frequency ujq, shifted by the cavity pulls xjp 
[3(1 [36| . Finally, the third term describes the effective 
qubit-qubit coupling, also called J-coupling or transverse 
exchange coupling, 

3 \ j 3 / 

a flip-flop interaction mediated by virtual photon ex- 
change. 

The transverse exchange coupling in Eq. ^ leads to 
an avoided level crossing of the excited qubit states fl9| . 

At qubit-resonance, where 6 q = u) g e — tu g 2 J = 0, the 

size of the splitting is 2J = h^2j ^Qj^d^/^j- The new 



cigenstatcs are the symmetric triplet states \gg), |ee) and 
\tjj s ) = (\ge) + \eg))/y/2, as well as the anti-symmetric 
singlet state \ip a ) = (\ge) — \eg))/y/2, see Fig. j2ja). In 
the maximally entangled states \ip s / a ) a single excitation 
is shared between the two qubits. More generally, for 
S q 7^ the eigenstatcs of the Hamiltonian in Eq. (TTJ) can 
be parametrized as 

\tp s ) = sm6 n \ge) + cos0 n |eg), (3) 
\ip a ) = cos9 n \ge) - sm9 n \eg), 

with the mixing angle 6 n determined by cos 29 n = 
-8 q /^AJ 2 + <5 2 and sin20„ = 2J/^AJ 2 + 5^. The sep- 
arable qubit states \eg) and \ge) are asymptotically real- 
ized, \ip a ) — > \eg) and \ip s ) — > \ge), for large qubit-qubit 
detunings [S q — > 00), as indicated in Fig. [^b). 

We have performed two sets of experiments using sam- 
ples with different parameters listed in table HI In these 
experiments, the energy spectrum of the coupled qubits is 
probed by monitoring the transmission through the res- 
onator while applying a second spectroscopy tone [17| at 
frequency w<j. For the spectroscopy measurement shown 
in Fig. [2jb) , the first qubit is kept at a fixed frequency 
ujg]) j 2tt and the second qubit frequency u> g 2 J / 2tt is swept 
across the avoided crossing by changing its flux bias using 
external coils. The value of J can be extracted from a fit 
of the upper and lower branch of the avoided crossing to 
the function 

/("; o>$,J) = (V + 4e) ± V (42 - + 4J2 ) /2. 

(4) 

where to is in this parameter regime an approximately 
linear function of the flux $ threading the second qubit 
loop. The fit parameters are the transition frequency 
ujg 1 } of the first qubit and the coupling strength J. Both 
are determined with a precision of typically better than 
0.5 MHz. In this particular example we find u> g e /2tt = 
5.210 ± 0.00005 GHz and J/2tt = 10.06 ± 0.06 MHz. 

Two additional features are observed in Fig. |5Jb). 
First, a third spectroscopic line centered between the up- 
per and the lower branch appears at higher drive powers. 
This is a signature of a two-photon transition from the 
ground state (\gg)) to the doubly excited state (\ee)) of 
the coupled qubit system that is only allowed directly at 
the anti-crossing. This is discussed in Section |VlJ Sec- 
ond, the upper branch shows a transition to a dark res- 
onance at the avoided crossing, which can be explained 
by the symmetry of the states with respect to the spec- 
troscopic drive, see Section IVl 

III. COUPLING TO THE FUNDAMENTAL 
RESONATOR MODE 

According to Eq. (0) the coupling J = %W g( 2 ' / A 
scales inversely with the detuning A = A' 1 ' = A' 2 ', con- 
sidering only a single resonator mode. We have recorded 
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FIG. 2: (a) Energy level diagram of two transversely coupled transmon qubits. (b) Spectroscopic measurement of the avoided 
level crossing in sample A as function of normalized flux 3>/<l?o threading the first qubit loop with the second qubit at a fixed 
frequency. The solid lines indicate energy levels calculated from a diagonalization the two-qubit Jaynes-Cummings Hamiltonian. 
(c) Experimentally extracted value of the coupling strength J as a function of qubit frequency (dots). Lines indicate calculated 
values of J(oj ge ) for different models, see text for details. 



the avoided crossing between the two qubits at different 
dctunings A from the fundamental mode of the resonator 
using a spectroscopic measurement performed on sample 
A. The corresponding parameters are listed in Table HI 
The measured values of J shown in Fig. |^c) are deter- 
mined for each detuning from a fit as described in Section 

im 

Considering only one relevant resonator mode and con- 
stant g, the strength of the inter-qubit coupling is ex- 
pected to be symmetric about the resonator frequency 
(Eq. (|2|); thin gray line in Fig. Hie)). The asymmetry 
in the data can partly be accounted for by including the 
frequency-dependence of the coupling g, as explained in 
Appendix [X] It follows from the transition matrix el- 
ements that J scales proportional to the transition fre- 
quency u! ge of the qubits. This scaling factor leads to 
an asymmetry of the exchange coupling around the res- 



Parameter 


Sample A 


Sample B 


CJO /2-7T 

k/2-k 

Ef/h 

E ( c ] /h 

Ef/h 

Ef/h 

gf/2-K 
gf/^ 


6.44 GHz 
1.57 MHz 

232 MHz 

233 MHz 
35 GHz 
38 GHz 

133 MHz 

134 MHz 


3.34 GHz 
1.91 MHz 
148 MHz 
153 MHz 
409 GHz 
375 GHz 
43 MHz 
42 MHz 



TABLE I: Parameters of samples A and B as determined from 
independent measurements, luo denotes the fundamental fre- 
quency, k the cavity decay rate, Ef 2 ^ the charging energy, 



-Bj 1 ' 2 ' the maximum Josephson energy and g^'*' the coupling 
strength to the fundamental cavity mode of qubits 1 and 2. 



(1,2) 



onancc frequency luq that improves the agreement with 
the data (dashed red line, Fig^c)). To check, whether 
the remaining discrepancy originates from the dispersive 
approximation, we have also done a numerical diagonal- 
ization of the full generalized Jaynes-Cummings Hamil- 
tonian (not shown) . This calculation agrees with the dis- 
persive model within the errors of the experimentally de- 
termined values of J. 

For a quantitative agreement, higher harmonics of the 
resonator have to be considered. The particular imple- 
mentation of the resonator as an open-ended coplanar 
waveguide supports higher harmonics at integer multiples 
of the fundamental frequency 34j, see Fig. EJa). Each of 
these higher modes provides a channel for the exchange 
of virtual photons between the qubits determined by the 
detuning and the coupling to the harmonic mode 
j as indicated in Fig.[3][b). Above the fundamental mode, 
the coupling to the first harmonic mode j = 1 contributes 
significantly to the qubit-qubit coupling, which results in 
an asymmetry with respect to the detuning A. Including 
four modes in Eq. © to determine the expected value 
of </, good agreement with data is obtained (thick green 
line, Fig. He)). 

It is important to also include the alternating sign 
of the electric fields at the qubits' position in the cal- 
culations. As depicted schematically in Fig. Ela), the 
electric field of the fundamental and higher even modes 
(j=0, 2, 4,. . . ) have always opposite sign, i. e. a relative 
phase of 7r, at either end of the microwave cavity, whereas 
odd modes (j=l, 3, 5,. . . ) have equal sign. Thus, higher 
harmonics add with different signs to the effective cou- 
pling strength. A priori, the sum in Eq. © has to be ex- 
tended over all modes. The sum does, however, not con- 
verge, as discussed also in the context of Purcell-limitcd 
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FIG. 3: (a) Spatial mode structure of a coplanar waveguide 
resonator. Qubits are positioned at opposite ends of the A/2 
resonator with coupling of alternating sign, — (— 
to the j-th resonator mode at frequency ujj . (b) Energy level 
diagram and coupling scheme for two qubits with transition 
frequencies edge and u) g V coupled to a transmission line cavity 
with fundamental frequency cjo. Energy levels with a photon 
in the j-th resonator mode (gglj ) or a qubit excitation (geO or 
egO) are shown. The exchange interaction J depends on the 
detunings A* 1 ' and the coupling strengths g^ of both qubits 
i = 1, 2 to the resonator mode j. 

qubit decay rates in [37| • The proportionality of the cou- 
pling g^gf oc ujj = (j + l)u>Q (Eq. IA2|) together with 
the same proportionality of the detuning Aj oc (j + 1)cjo 
for large j leads to an non-converging series alternating 
between the two values 

J even = ft^(-1)^+ 1 ^ and J odd = n J2(-iy +1 ^-. 

3=0 3 j=0 3 

Apparently, a cut-off frequency has to be imposed to ob- 
tain physical results. In Fig.^c) we have also included a 
plot of Jeven when terminating the sum at the fourth har- 
monic (k = 2) (dotted blue line). It is observed that the 
difference between an even and odd number of modes 
is significant, up to 25% of the coupling strength. We 
have also verified that this is not an artifact of the dis- 
persive model by numerically diagonalizing the Jaynes- 
Cummings Hamiltonian. 

In our measurements we observe enhancement of the 
exchange coupling, inversely proportional to the detuning 
of the qubits to the cavity mode. The asymmetry around 
the mode is attributed to higher harmonic modes that 
contribute to the measured (renormalized) J. However, 
to compute the coupling strength the number of included 
modes has to be restricted by imposing a high-frequency 



cut-off. Physically, there are several mechanisms conceiv- 
able. The energy needed to overcome the pairing inter- 
action of Cooper pairs sets an upper frequency of about 
700 GHz for Niobium. The electric field across the trans- 
mon averages out when the wavelength of the photons 
becomes comparable to the size of the transmon at a fre- 
quency of about 400 GHz. Also, radiation or dielectric 
loss mechanisms and the photon loss rate through the 
coupling capacitors increase at higher frequencies [37| . 
Current experiments arc, however, not designed to work 
at frequencies higher than approx. 15 GHz. The explo- 
ration of the relevant frequency range will require an elab- 
orate circuit architecture and will be challenging with 
current technology. 



IV. MULTI MODE COUPLING 

To assess the coupling to higher-order modes, we mea- 
sure the qubit-qubit coupling strength J as a func- 
tion of qubit transition frequency in a second sample 
B over a broader frequency range. In this sample the 
frequency of the fundamental resonator mode is lower, 
u>o/2n = 3.34 GHz, and the maximal Josephson energy 
Ej is higher by about one order of magnitude (sec Table 
U- sample B). As a result, the qubit transition frequency 
can be swept over several resonator modes. 

The measured inter-qubit coupling strength shows an 
enhanced value around each harmonic mode, as well as an 
overall increase with frequency (Fig.0]) in agreement with 
the discussion in Sec. lIIII A calculation based on the dis- 
persive Jaynes-Cummings Hamiltonian can explain the 
data qualitatively. However, including an even (N = 6, 
solid red line) or an odd (N = 7, dotted blue line) number 
of resonator modes yields again significant differences in 
the calculated value of J, neither of the curves resulting 
in good quantitative agreement. 

The resonator modes discussed so far cannot fully ex- 
plain the measured inter-qubit coupling. Large devi- 
ations located asymmetrically around ~ 8.5 GHz and 
« 12 GHz (Fig. [3} hint at the presence of an additional 
coupling mechanism at these intermediate frequencies. 
The measurement reveals that the coupling strength is 
asymmetrically modified between every two resonances 
alluding to an anti-resonance that mediates a qubit-qubit 
coupling channel of similar magnitude as the coplanar 
waveguide resonance. This can lead to an enhancement 
or suppression of the qubit-qubit coupling at intermedi- 
ate frequencies due to interference between multiple res- 
onances. 

To account for these spurious resonances, additional 
field modes are included in the model. Physically, these 
modes may be identified as the slotline modes of the 
coplanar waveguide, a differential excitation of the left 
and right ground plane [38j . Air-bridges or wire-bonds 
that connect the ground planes could effectively sup- 
press these modes, but have not been implemented in 
this sample. Treating these modes cquivalcntly to the 
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FIG. 4: Exchange coupling J versus qubit frequency in sample B. The vertical lines indicate the frequency of the coplanar 
waveguide modes ujj. Experimental data (dots) confirm the expected increase of J with increasing frequency. The solid red 
(dotted blue) line indicates the calculated J including N = 6 (TV = 7) resonator modes. The dashed green line is a fit to a 
model with additional resonances. 



coplanar waveguide modes in the derivation of the dis- 
persive multi-mode Jaynes-Cummings Hamiltonian ([1}, 
an extra contribution to the exchange coupling emerges, 



Ju 



J+J = J 



1 M 

1 ~(1)~(2) 



A 



(i) 



A 



(2) 



Jtot is then fitted to the data in Fig. @] assuming that 
the coupling strengths of both qubits are equal in magni- 
tude, | = \g\ i |, but have - like the coplanar waveg- 

u(i)i - (-i)l+i\gf)\ W ith 



uide modes - alternating sign, \g t 

I = 0, 1, For the fit we take four extra modes into ac 

count (Z = 1,2,3, 4) from which we obtain the resonance 
frequencies £1,2,3,4/2^ w {5.2, 8.4, 11.9, 14.8} GHz. 
Note, that these additional modes are not observed in 
simple transmission measurements [34| of the resonator 
and their frequencies can therefore not be determined 
independently. In the same fit the ratios gj/gj are deter- 
mined to £1,2,3,4/01,2,3,4 = {1(1), 0.5(2), 0.8(2), 0.7(1)}, 
which shows that the coupling strengths to the spurious 
modes gi and to the coplanar waveguide modes gi are sim- 
ilar in strength. This similarity hints at a highly localized 
field of the spurious mode with small effective mode vol- 
ume. Also, the relative sign between the couplings g( 

~(2) 

and g\ of the qubits to the spurious mode alternate with 
the mode number I. This implies that the field of the spu- 
rious modes has also either equal or opposite direction at 
the position of the qubits, like the coplanar waveguide 
mode outlined in Fig. [3] The qualitative agreement to 
the measured values of J is considerably improved by the 
inclusion of these extra modes, see dashed green line in 
Fig.|U 

The measured values of the exchange coupling J 
demonstrate the sensitivity of the qubit-qubit coupling 
to the full mode structure of the circuit. While a single- 
mode model is sufficient around a single resonance, quan- 
titative predictions require complete knowledge of de- 
signed and spurious resonances. Appropriate circuit de- 
sign and use of wire-bond or air-bridge connections of 
ground-planes on the chip can short out spurious modes. 



In contrast, additional resonances can also be incorpo- 
rated on purpose into the circuit design ,39| to modify 
the qubit-qubit coupling at certain frequencies. 



V. DARK STATE 

A characteristic feature of the avoided level crossing 
is the observation of a dark resonance, where the transi- 
tion from the ground state to the upper energy branch is 
forbidden and no signal is observed in spectroscopy mea- 
surements, see Fig. [5Jb). In fact, the symmetry of the 
states at the avoided crossing leads to a selection rule 
with respect to the spectroscopic drive [III l3(i| 



d = e(gWaW/AU+gW<r?/AW 



h.c. 



(5) 



through the resonator. To see this, we decompose the 
eigenvalues of the dispersive Hamiltonian ([1} into triplet 
states and a singlet state. These are the eigenstates of the 
permutation operator (Eq. lBl|) to the eigenvalue ±1. Ex- 
plicitly, the triplet states \gg) (ground state), |ee) (doubly 
excited state) and the symmetric state (\ge) + \eg))/\/2 
span the symmetric subspace, whereas the singlet state 
(\ge) — \eg))/\/2 is anti-symmetric under permutation of 
the qubits. This is equivalent to a decomposition of the 
system into a spin-1 and a spin-0 particle. The sign of 
the coupling constants g^ determines the drive symme- 
try. Equal (opposite) sign of the electric field at the po- 
sition of the two qubits leads to the positive (negative) 
sign of the second term in Eq. ([5]) and, consequently, to 
a(n) (anti-)symmetric excitation. For equal sign of the 
couplings, the drive term and the permutation operator 
commute (symmetric drive). Then, only transitions be- 
tween states of same symmetry are allowed (see Appendix 
|B|) and the anti-symmetric state stays dark at zero de- 
tuning. Vice versa, for opposite sign of the couplings 
the now anti-symmetric drive can connect symmetric to 
anti-symmetric states. 

In our experiments the symmetry of the drive is de- 
termined by the frequency of the microwave signal and 
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the distance d between the qubits. The relative sign of 
the field at the qubit positions Xi is determined by the 
phase difference A<fi = 4>{x\) ~ 4>{ x 2) = ^dd/c c g of the 
travelling wave between the qubits. Here, we have used 
the dispersion relation kd = 0Jd/c e g with the propaga- 
tion velocity c e g of light in the transmission line. As the 
qubits are located at the end of the transmission line res- 
onator, d is approximately the length of the resonator and 
sign changes happen at frequencies ui s = sirc e s/(2d) with 
s = 1, 3, 5, . . ., in between two resonances as indicated by 
the dotted lines in Fig. O^a). 

Whether the eigenstate with lower or higher energy is 
dark, depends - in the simplest model with only a sin- 
gle dominant resonator mode - on the qubit-resonator 
detuning. The higher (|V'+)) an d lower (\ip-)) en- 
ergy eigenstate of the Hamiltonian (|TJ) at the avoided 
crossing can according to Eq. ([3]) be written as ip± = 
(\ge) ± sign(J)|eg)) /v2- Note, that in this notation the 
subscript denotes higher (+) or lower (— ) energy and 
not the symmetry of the state. In Fig. [5] the energy lev- 
els of the coupled resonator-qubits system arc plotted as 
a function of the qubit transition frequencies, which are 
kept equal, along with their respective symmetries. For a 
coupling to the first harmonic mode J = jf /Ai < 
below the mode (Ai < 0, g[^g^ > 0) and J > above 

the mode (Ai > 0, g^ > 0). Consequently, below 
the first harmonic the higher energy state ip+ is the anti- 
symmetric singlet state. This state cannot be excited 
from the ground state with a symmetric drive and stays 
dark (Fig. \5j[b) - sub-panel A) . At the avoided crossing 
with a resonator mode the lower and higher energy qubit 
state swap their symmetry due to the sign change of the 
detuning. Since the drive does not change its symmetry, 
also the dark and bright state energies are interchanged 
and the dark state appears at the lower branch (Fig.[5jb) 
- sub-panel B). The dark state is always closer in fre- 
quency to the resonator transition. 

If the second harmonic mode dominates the coupling, 
the situation is reversed since the coupling constants 
have different signs. Below this mode, the ip+ (higher 
energy) state is symmetric, and ip- (lower energy) is 
anti-symmetric. Still, the dark state appears at the 
upper branch of the avoided crossing (Fig. [5][b) - sub- 
panel C) and switches to the lower branch above the res- 
onator (Fig. [5jb) - sub-panel D). The reason is that the 
drive symmetry changes as well in between two resonator 
modes as explained above, implying that the drive field 
changes from a symmetric (around the first harmonic 
mode) to an anti-symmetric drive (around the second 
harmonic mode). The drive can then induce transitions 
between the ground and an anti-symmetric state, but not 
to the symmetric state i[>+ of the upper branch. The var- 
ious conditions leading to the identification of the dark 
state branch are summarized in table |TTJ In particular, if 
the drive (line 4) has different symmetry than the higher 
energy state ip+ (line 3), ip-y. remains dark. 
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TABLE II: Symmetry considerations leading to the dark state 
at the lower or upper energy branch of the avoided level cross- 
ing. © (©) denotes a positive (negative) value or symmetry. 




anti-symmetric state 
symmetric state 



3.2 6.4 6.6 6.8 7.0 9.6 9.8 10.0 10.2 10.4 

joint qubit transition frequency [GHz] 



(h) 




FIG. 5: (a) Energy levels of the qubits-resonator system when 
the qubits are degenerate and their transition frequencies are 
simultaneously swept across the resonator modes. The param- 
eters are those of sample B (Table H]). (b) Spectroscopic mea- 
surements of the anti-crossing in the regions labeled A, B, C 
and D in (a). 



VI. TWO-PHOTON TRANSITION 

The spectroscopic line between upper and lower branch 
of the avoided level-crossing in Fig. [UJb) is a two-photon 
transition from the ground state \gg) to the doubly ex- 
cited state |ee). Similarly, this transition has also been 
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observed in a phase qubit coupled coherently to a two 
level fluctuator in the tunnel barrier of the Joscphson 
junction comprising the qubit and in molecular spec- 
troscopy of two nearby molecules [4lJ . It becomes visible 
only at the center of the avoided level crossing, again a 
manifestation of the symmetry properties of the system. 
The rate of the corresponding two-photon transition [42| 



r= 2n 



(ee\H d \m)(m\H d \gg) 

^ 0J m -0J d 
m 



is given by a sum over the intermediate states m. Off 
the avoided crossing the qubits are effectively decoupled. 
In this case the intermediate states are m = ge, eg. The 
transition is then prohibited due to destructive interfer- 
ence between the two possible paths, gg <-> eg O ee or 
gg <-> ge -H> ee, connecting the ground to the doubly ex- 
cited state. Due to the opposite sign of the detunings 
Wge/eg — u d , the two terms in the sum cancel and the 
transition rate T = 0. With the qubits at resonance, 
the intermediate states are m = ip + ,tp~ and one term 
in the sum vanishes due to the forbidden transition to 
the dark state. With only one possible path connecting 
the gg to the ee state, no interference takes place and 
the transition becomes allowed. The enhanced transition 
rate can be employed for directly creating the maximally 
entangled state \<f>) = (\gg) + \ee))/y/2. 



VII. CONCLUSION 



formed by the microwave resonator, dissipativc dynamics 
of separable and entangled states can be studied. The ex- 
change coupling can also be useful for building two-qubit 
gates when fast flux-pulses are applied to tune the qubits 
into resonance. In the context of quantum information 
processing, the resulting SWAP gate forms a universal 
two-qubit gate with short operation times. Moreover, 
the transverse exchange coupling mechanism described 
in this article mediates interaction not only between two, 
but an arbitrary number of distant qubits, an interest- 
ing playground for studies of collective phenomena with 
superconducting circuits, where the interaction is not re- 
stricted to nearest-neighbours. 
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Appendix A: Frequency-dependence of the 
qubit-resonator coupling 

The coupling g to the resonator is proportional to the 
rms voltage fluctuations of the vacuum field V r ° ms at the 
position of the i-th qubit and to the off-diagonal matrix 
element [HI H3] of the charge operator h, 



g = 2e/31/° ms (ff|n|e). 



(Al) 



We have analyzed the coupling between two distant 
qubits mediated by the harmonic modes of a resonator. 
We have observed an overall increase of the exchange 
coupling with frequency as expected from a model in- 
cluding higher-harmonic modes of the coplanar waveg- 
uide resonator. Good qualitative agreement over a wide 
frequency range between the dispersive model and experi- 
mental data is obtained when taking spurious resonances 
of the coplanar waveguide in addition to the coplanar 
waveguide modes into account. Hence, measurements of 
the transverse inter-qubit coupling can be employed to 
detect and investigate spurious global coupling channels 
between distant qubits, complementary to measurements 
of single qubit spectra used to detect spurious local res- 
onances (43l |44| . How many higher harmonic modes to 
include in the theory, i. e. where to set a high-frequency 
cut-off, can, however, not be decided on the basis of cur- 
rent measurements. 

In addition, we have observed dark states and en- 
hanced two-photon absorption at the avoided level cross- 
ing in spectroscopic measurements. These characteristic 
features are based on the relation between the symmetry 
of the drive and the singlet and triplet states formed by 
the coupled qubits, which also explains the dark state 
at either the lower or higher energy branch. These sym- 
metries also affect decay processes of singlet and triplet 
states and, together with the non-trivial environment 



The prcfactor f3 is determined by the geometry of the 
circuit used in our experiments. In the large Ej/Eq 
limit, realized in the devices, the matrix element is pro- 
portional to the square-root of the qubit transition fre- 
quency, (<?|".|e) oc y/ojge- The vacuum field V^ ns is 
proportional to the s quar e root of the mode frequency, 

= + i)wo BM- 

For a single qubit on resonance with the j-th res- 
onator mode, uige = (J + 1)cjq, the scaling of the qubit- 
resonator coupling is approximately linear in the mode 
number, g oc (j + 1)(Jq. To verify the linearity of the 
coupling strength we have measured the vacuum Rabi 
splitting of a single qubit up to the third harmonic res- 
onator mode (Fig. HJa)) using qubit 2 of sample B (for 
the parameters, see Table [I]). The simple estimate shows 
good agreement with the measured coupling strengths 

ffo 2 i,2 a/ 271 " = i 42 : 84 > 125 > 162 l MHz (dashed line in 
Fig-lSla)). The parameter (3 = 0.20, obtained from a lin- 
ear fit to the analytic model in the large E j / Eq limit [33| , 
agrees with the designed value within 10%. A numerical 
simulation of the transmon including four energy levels 
can explain also the slight deviations from the linear de- 
pendence at high frequencies (solid red line, Fig. IHa)). 

In the case of two resonant qubits with identical cou- 
pling strength to the j-th harmonic mode, the exchange 
coupling strength is, according to Eq. (|2|), linear in the 
qubit transition frequency uj ge and the mode frequency 
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1 2 
mode number i 



(b) 



-0.2 



99 

ee 



5 q [GHz] 
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H d P + PH, 



(B2) 



This relation can be fulfilled only if the drive transforms 
a symmetric (tp s ) to an anti-symmetric (ipa) state, or vice 
versa, such that 



FIG. 6: (a) Coupling strength gj 2 ' of qubit 2 to the j-th 
harmonic mode in sample B. (b) Expectation value of the 
permutation operator P for the eigenstates of the Jaynes- 
Cummings Hamiltonian as a function of qubit-qubit detuning 

S„. 



[H d ,P] + ^ s 



H d Pip s + PH d ip s 

P d l\) s + Plp a 
4>a-4>a= 0. 



(B3) 



(j + l)wo, 

J <x 9j9j ~ 9j « u ge Uj = {j + l)uj ge u . (A2) 

When the detuning of the qubits to the resonator mode is 
varied, the frequency (j + l)wo of the j-th resonator mode 
is constant and J scales proportional to the transition 
frequency uj ge . 

Appendix B: Dark state symmetry 

The spectroscopic drive H d oc (a^ — <r^) + h.c. [3(3| 
anti-commutes with the permutation operator 

P = ( CT <V 2) + 4V 1 )) + (l + aPaW) /2, (Bl) 



The symmetry of the coupled qubit states (see Figure 
[2^a)) can be characterized by the corresponding expecta- 
tion value of the permutation operator (P) . (P) is one for 
the \gg) and the |ee)-state, i. e. the \gg) and |ee) states 
are symmetric for all dctunings. For the symmetric and 
anti-symmetric states ip s and ip a formed at zero qubit- 
qubit detuning, (P) is 1 or — 1 indicating that these states 
are eigenstates of P with well-defined symmetry. For 
non-zero detuning between the qubits, S q ^ 0, the eigen- 
states of the Hamiltonian (JTJ) do not have well-defined 
symmetry and (P) approaches asymptotically zero for 
large dctunings (Fig. El(b)). Hence, no strict selection 
rules are imposed off the level crossing and the transition 
between ground state and single excited states is allowed. 
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